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Abstract 



We investigate the potential effects of the scalar flavour changing neutral currents that 
are generated e.g. in supersymmetry with tan /9 ^ 1 in the b sl~^l~ transitions. Using 
the experimental upper limit on BR{B^ — ^ f^^f^~) we place stringent model independent 
constraints on the impact these currents may have on the rates BR{B — > Xg/i+yU") and 
BR{B Kfi^fi~). We find that in the first case, contrary to the claim made recently in 
the literature, the maximal potential effects are always smaller than the uncertainty of the 
Standard Model NNLO prediction, that is of order 5-15%. In the second case, the effects 
can be large but the experimental errors combined with the unsettled problems associated 
with the relevant formfactors do not allow for any firm conclusion about the detectability 
of a new physics signal in this process. In supersymmetry the effects of the scalar flavour 
changing neutral currents are further constrained by the experimental lower limit on the 
B^-B^ mass difference, so that most likely no detectable signal of the supersymmetry 
generated scalar flavour changing neutral currents in processes B Xsfi^fi~ and B — >■ 
Kfi'^fi~ is possible. 



1 Introduction 



Rare processes involving the 6-quark, intensively studied at present in several experiments 
(BaBar, BELLE, Tevatron), play an important role in supersymmetry (SUSY) searches 
via virtual effects of the new particles. This is because in the minimal supersymmetric 
extension (MSSM) of the Standard Model (SM) the Yukawa couplings of the 6-quark 
to some of the superpartners of the known particles and/or to the Higgs bosons can be 
strong enough to produce measurable effects. A celebrated example is the radiative decay 
B Xsj whose experimentally measured rate P agrees very well with the SM prediction 
|2] and, consequently, puts constraints on the MSSM parameter space. These constraints 
become particularly stringent if the ratio Vu/vd = tan (3 of the vacuum expectation values 
of the two Higgs doublets is large, that is when the coupling of the right-chiral 6-quark to 
charginos and the top squarks is enhanced: agreement with the experimental value can 
be then obtained either if all these sparticles as well as the charged Higgs boson are 
sufficiently heavy (in which case there is little hope to detect their virtual effects also in 
other rare processes), or if the virtual chargino-stop contribution to the 6 ^ 57 amplitude 
cancels against the top-charged Higgs boson contribution. The latter solution requires of 
course a certain amount of fine tuning, which becomes, however, of tolerable magnitude 
for Mh+ > 200 GeV and sparticles weighting not less than a few hundreds GeV. 

A very interesting feature of the large tan/3 SUSY scenario is the generation at one 
loop of the (tan^ /5)-enhanced flavour violating (FV) couplings of the neutral Higgs bosons, 
AP (the CP-odd one) and (the heavier CP-even one), to the down- type quarks [3]. 
Being operators of dimension four, these couplings remain unsuppressed even for heavy 
superpartners of the known particles (gluinos, squarks, charginos). If the flavour violation 
is minimal (the so-called MFV SUSY), that is if the Cabbibo-Kobayashi-Maskawa (CKM) 
matrix is the only source of flavour and CP violation, the FV couplings of and are 
very sensitive to the mixing of the left and right top squarks. (Induced by these couplings 
FV decays of the neutral MSSM Higgs bosons have been investigated in ref. jl].) The 
exchanges of the neutral Higgs bosons generate then | AF| = 1 [30111 and | AF| = 2 [Slini 
dimension six operators which contribute to the h sl^l~ and hs hs \^ transitions. 
For and not much heavier than the electroweak scale these operators, called because 
of their Lorentz structure the scalar operators, can significantly change the predictions of 
the SM. 

Phenomenological consequences of the scalar operators have been analyzed in several 
papers 0^11311011111113121113113113113111111312011^1^ both in supersymmetry 
with minimal (MFV) and nonminimal flavour violation. In particular, it has been shown 
|3II3II3IIlllini that even in the MFV SUSY the effects of the scalar operators originating 
from the FV couplings of and A^ can, for large mixing of the left and right chiral top 
squarks, increase BR{B^ — > fi'^fi~) and BR{B^ /^^yU~) by 3 — 4 orders of magnitude. 
The upper bound on the first of these branchin fractions set recently by CDF j2Sl 

BR{B°^ /i+/i") < 0.95 X lO"*' at 90% C.L., (1) 

(which improves the previous limit BR{B^ /U^/W ) < 2 x 10~^ [24j) puts therefore on the 
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MSSM parameter space a nontrivial constraint which is to a large extent complementary 
to the one imposed by the measurement of BR{B X^j). On the other hand, as 
shown in a measurement of the fJ'^f^' signal at the Tevatron Run II, possible 

if BR{B^ — >■ fJ.~^fi^) ^ 2 X 10^®, would rule out such models of the soft SUSY breaking 
terms generation like anomaly and gaugino mediation as well as gauge mediation scenarios 
with low messenger scale and small number of messenger fields. 

The impact of the FV couplings of and on the |AF| = 2 transitions 5^ 
5° ^ 50 ^as analyzed within the MFV SUSY in refs. HZl CHI 1201 ■ It was found that 
the contribution of the |AF| = 2 scalar operators constructed out of these couplings to 
the amplitude of the B^-B^ mixing is negative and can be very large {B^-B^ mixing is 
affected negligibly). Part of the parameter space corresponding to tan /? 3> 1, light and 

and substantial stop mixing, allowed by the experimental limit on BR{Bg fJ'^f^~) 
then available, was eliminated by the condition that the calculated B^-B^ mass difference 
AMs is not smaller than the experimental lower bound AMs ^ 14/ps fl^. Even with 
the new bound ((T)) the constraints on the MSSM parameter space imposed by the B^-B^ 
mixing are in some cases stronger than the ones stemming from the dimuon channel. 

The effects of the scalar operators in the exclusive transitions B —* K^^^' and 
B —>■ K*fi~^fi~ have been investigated in ^OJEj. Their impact on BR{B —>■ K*fi~^fi~) 
has been found to be very small. On the other hand, potential effects of the scalar 
operators in — > Kfi'^fi" could be quite sizeable in principle, but the experimental limit 
BR{B^ K~^fi~^fi~) < 5.2 X 10~^ available at that time was too weak to provide 
constraints stronger than the experimental upper limit for BR{B^ — > fi~^fi~). Finally, 
the effects of the scalar operators in the inclusive decay rate BR{B — >■ Xsii^^~) have 
been taken into account in several papers devoted to general investigation of the potential 
SUSY effects in radiative B decays or in the studies of the specific SUSY scenarios like 
the minimal SUGRA, but have not been directly confronted with the bounds provided by 
the B^ decay and B^-B^ mixing. 

In this paper we fill this gap. We begin in section 2 by recalling the NNLO predictions 
of the SM for BR{Bg — > and BR{B Xsfi~^^~) improving slightly in the latter 

case the estimates of the theoretical uncertainties compared to those given in ref. pSj . 
Then in section 3, following ref. ^^1; we asses in a model independent way how big effects 
of the scalar operators in the BR{B —>■ Xsfi'^fi") and in BR{B —>■ Kfi^fi^) decays are 
still allowed by the CDF bound BR{B^ f^'^fJ''') < 0.95 x 10"^. We show in particular, 
that the huge effects of the scalar operators found recently in BR{B —>■ XsH'^n~) in 
ref. 1221 excluded by these constraints. The results of section 3 are valid generally, 
independently of the mechanism that generates the scalar operators. Finally, in section 4 
we concentrate on scalar operators in the MFV version of the MSSM (in which the squark 
mass matrices are aligned with the quark ones - see [201 for more detailed explanations) 
and specify the maximal effects of the scalar operators in BR{B Xsfi~^fi~) and in 
BR{B — > i^/i^/i") allowed by the experimental limits on both, the B^ fi~^fi~ rate and 
the B^-B^ mass difference. We summarize the situation in the last section. 



2 



2 b ^ sl^l and b dl^l transitions in the SM 



Under the assumption of minimal flavour violation, the effective Hamiltonian describing 
the h — > sl^l~ {b —>■ dl^l~) and b —>■ sj transitions takes the form jSU] 

n''^ = -2V2G^V,f*V,f{f:Cx{^^)Ox{f^)+ E E (2) 

\X=1 l=e,ti,T x=s,p J 

with the following set of operators Ox 1^ [T3] 

0,, = {sl1^T^cl){cl1,T%l) 

Os = isLYbL) E (W) 
q=u,d,s,c,b 

q=u,d,s,c,b 

q=u,d,s,c,b 

q=u,d,s,c,b 

O7 = ^{sLa^''bji)F,, (3) 
9; 

Os = -irsLa^''T%n)G% 
9s 

9s I 

Ol0 = ^{sLYbL)J2{lY7'l) 

9s I 
9s 

0'p=^isLbR){h'l) 

9; 



and Oiu, O2U obtained from Oic and O2C by the replacement c u, and the Wilson 
coefficients Cx{lj) organized as [301 

Cx{,) = cPi,) + '^2Ci\f^) + '^.cPi,) + ... (4) 

The coefficients Cx computed at some scale /xq ~ are subsequently evolved down to the 
scale fib ~ rrib, where their matrix elements between the hadronic initial and final states of 
the process under investigation are computed either by lattice methods or perturbatively 
to the required accuracy in as{^b) = 9s{f^b)/4:7i. At the matching scale /xq only the 
coefficients of the operator O2 starts at order (q;s)°; for the remaining ones Cx\l^o) = 0. 
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2.1 11+ ij- in the SM 

In the SM the Wilson coefficients Cs and Cp are neghgible and the only operator relevant 
for the 5° ^ — > transitions is Oiq. Its Wilson coefficients Ciq' and C^g'* at the matching 
scale are known 1^- Since the quark part of Oio is a (partially) conserved chiral 
current, the QCD evolution of Cio is simple, i.e. Cio(yUfe) = [Q;s(yUb)/as(/io)]C'io(/io)- This 
leads to the well known prediction 
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J' 
TT 



47rsin 9w I 



2 

TT) 

I - ^j^\v:M\nxt)\\ (5) 



where 



-i^Yixt) = C[l\x,) + ^f^Cg)(x„/io) (6) 
sm Ow Idtt^ 

and Xt = (m^^(/io)/Mvi/)^. C[Q'{xt) is given by the function Yo{xt), which can be found 
e.g. in and CiQ'{xt, Ho) has been computed in [221 (it can be also extracted from [HT]). 
For mp(mt) = (166 ± 5) GeV, a,(M^) = 0.119 and using /iq = = 174.3 GeV 

Y{xt) = r] (0.971 ± 0.046) (7) 

where r] = 1.01 accounts for the effects of C[q' . For sin^ 6w = 0.23124 and Oem = 1/128 
this gives 



i?i?(i?° ^ /xV") = (3.64 ± 0.33) X 10-^ X (j^^) 
5i?(5^ ^ /X+/X') = (1.39 ± 0.13) X 10"^° y ' > 




1.542 ps/ V 203 MeV/ \ 0.009 



of the SM predictions (of order ~t24% ~t3o% in the case of the 5° and decays. 



where the errors correspond to the variation of m^^(mt). The dominant uncertainities 



respectively) come from the factors Fb, = (238 ±31) MeV and Fb, = (203 ± 27i° 



20) 



MeV jSl] that parametrize the nonperturbative hadronic matrix element of the Oio op- 
erator. The uncertainty associated with Amf^^{mt) = 5 GeV, with the electromagnetic 
corrections and, in the case the B^ decay with the value of \Vtd\, are much smaller. 

The corresponding branching ratios for the e+e~ channel are suppressed by the factor 
(mg/m^)^ ~ 2 X 10~^ and, hence, unmeasurably small; those for the t~^t~ channel are 
enhanced by (m^-Zm^)^ ~ 283 but taons are very difficult to identify experimentally. 

The present experimental bounds BR{B'^ l^+l^') < 0.95 x 10"^ ^ and BR{B^ 
fi~^fj'~) < 1.6 X 10~^ [331122] are 3 orders of magnitude above the predictions ^ and still 
leave a lot of room for new physics. 
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2.2 The inclusive process B XJ^l in the SM 



The general formula for the differential width of the B — > Xgl^l decay reads I29j: 



ds 



7687r5 ' 



V:M^\''\l,r,,s)X'l\l,rJs,rJs) 



+h{s)G2{s) \cf{s,ii,)'^ + /4(s)G3(s)Re (cf (s, /i,; 



(9) 



where s = /ml is the "reduced" invariant mass of the lepton pair and 

A(a, 6, c) = + 6^ + - 2a6 - 2ac - 2bc . 



.(1)*/ 



(10) 



The function Gc{s, Ai, A2) accounting for the l/m"^ nonperturbative contribution has been 
found in The l/mf nonperturbative contributions summarized by the functions 

Gi{s, Ai, A2) have been calculated using the heavy quark expansion technique in P7| 1^. 
The functions Gc{s, Ai, A2) and Gi{s, Ai, A2), which depend on the parameters Ai ~ —0.2 
GeV^, A2 = 0.12 GeV^ are given in eqs. (29-31) of W^. Finally/ 

fi{s) = s{l + r, - s)A(l, rs/s, rjs) + (1 - r, + s)(l - r, - s){l + 2rjs) 
+ 6n(l + rs- s) 

/2(s) = s(l + r, - s)A(l, rjs, rjs) + (1 - + s)(l - - s)(l + 2r,/s) 
- 6n(l + - s) 



/3(^) 
/4(^) 

/5(^) 

/6(^) 
/7(^) 



(4/s)(l + 2r,/2) [2(1 + r,)(l - r,f - s{l + 14r, + 

I2{l + 2rjs) [{l-Tsf -s{l + r,) 
3 



(11) 



1 + - s)(s - 4n 



'1 + r. - s)s 



6^/^7(1 - - s) 



where r/ = mf/ml, Vg = ml/ml. In the NNLO approximation the coefficients CY^{s,fj.b), 
Gg^{s,fib) and ClQ{s,fib) summarizing the effects of the QCD running from the scale 
fiQ ~ rrit down to the scale /if, ~ rrib and the matrix elements of the relevant operators 
from the list © can be compactly written as 

Cf{s,^^b) = (l + ^u;,{s)\ A, 



^The functions fz{s) and /7(s) differ from the corresponding expressions in ref. (29) . Due to the extra 
piece —s^ the function /3(s) as given here reproduces in the limit rris = the result obtained in earlier 
papers for the coefficient of |C'|'^p. We also confirm that the sign of /7(s) is as in the earlier papers jJS] 
(opposite to the one in UHl)- 
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Cf{s, fi,) = (l + ^^u,{s)] \A, + ng{ml/ml s) + U,g{l, s) + W,g{0, s)] (12) 



TT 



4:71 



(crFr(.)+c7fFf(.)+c«Ff(s 



as(/ife) 



TT 
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where Ai, Tg, Ug, Wg, the function g{z, s) can be found in and the exphcit formulae 
for the functions Fj^\s) and uji{s), vahd for s ^ 0.25 are given in refs. [SHl-^ Wilson 

coefficients C[^\ 6*2°^ ^^id C^^^ can be found e.g. in eqs. (E.9) of (41 j. One should also 
remember to expand the formula Q only up to terms of order ^^(/ib) and to replace 
ujt^s) and ujg{s) by uj7g{s) in the interference term. Inclusion to C^^, Cg^ and Cfg of 
the 0(as(/ib)) corrections to the matrix elements^ of the relevant operators significantly 
decreases the dependence of the final result on the renormalization scale fib [3SI- Since 
similar 0{as{fib)) corrections to the matrix elements of the operators O^^ p are not known 
at present, their contribution to the rates of the inclusive B — >■ Xsl^l~ processes has the 
uncertainty associated with choice of the scale /i;, larger than do the contributions of the 
remaining operators. 

In oder to get rid of the factor in the formula Q not introducing at the same time 
large uncertainty associated with the value of the charm quark mass we follow the trick 
proposed in ^] and normalize the rate to width of the charmless semileptonic decay 

dBR{B X,/i+/i-) BR{B X^eu^) 



ds C 
where the function h{z) is given e.g. by the formula (48) of [HI] and the factor C 

Vub 



c = 



cb 



TjB ^ X.ez/e) 



has been calculated in gl]: C = 0.575 x (1 ± O.Olport ± 0.02ai ± 0.02a) = 0.575 x (1 ± 
0.03). To remain conservative we will double this uncertainty and use C = 0.575 x 
(1 ± 0.06). The poorly known nonperturbative parameter Ai approximately cancels out 
between the numerator and the denominator. With this trick the residual dependence 
on z = ml/ml is negligible for ml/ml varying between 0.27 and 0.31; the uncertainty 

^Complete results for the matrix elements, valid in the entire range of s, have been reported in pHj 
but are not yet publicly available. 

^In this analysis we neglect the contribution of the real gluon bremsstrahlung calculated in which 
changes the result by ~ 1%. 



6 



of the differential branching fraction arising from the normahzation is then dominated 
by the ~ ±6% uncertainty of the factor C fll4|) . It is therefore much smaller than the 
uncertainty of order ±15% attributed to the differential branching fraction normalized 
directly to BR{B XcCVe) in ref. jlHj by varying ml/ ml in the range 0.25—0.33. 

The dominant source of uncertainty remains the dependence on yU;, which for s < 0.25 
is estimated (by changing /i^ between 2.5 GeV and 10 GeV) to be of order ±7% [39^. Of 
comparable magnitude can be however also the uncertainty related to the electromagnetic 
corrections to the running (and their mixing with others) of the Og and Oiq operators, 
which is unknown at present."^ Simple estimate of this effect is obtained by varying Oem in 
the formula ^ between 1/128 and 1/133. This suggests additional ~ 8% uncertainty of 
the predicted branching ratio. Finally, the parametric uncertainty related to the variation 
of mp(mt) = (166 ± 5) GeV is of order ±(6 - 7)%. 

The differential rate Q can be integrated over various domains of s. The most reliable 
theoretical predictions are obtained for 0.05 < s < 0.25 because for this range the nonper- 
turbative effects associated with the cc resonances are small and the NNLO calculation is 
complete. For this region, using m^^(mt) = 166 GeV, m^, = 4.8 GeV, aem = 1/128 and 
\VtsVtl/V,b\ = 0.976 we get: 

BR{B ^ X,/iV")o.o5<s<o.25 = (1-46 ± 0.11 ± 0.10) x lO"'' (15) 

where we have used BR{B — ^ XceUf,) = 0.102. The first uncertainty comes from the 
fib dependence and the second one from Amf^^{mt) = 5 GeV. To this one has to add 
the 6% uncertainty from the C factor and (conservatively) a ~ 8% uncertainty from 
the electromagnetic corrections. Adding all these uncertainties in quadratures we finally 
assign to the result the uncertainty of order ±14%. 

Integrating the differential rate (jHI) over the entire domain^ Smin < s < Smax where 
■Smin = 4mf/m^, Smax = (1 ^ rns/mbY one obtains the so-called "nonresonant" branching 
fraction which can be compared with the experimental data provided the contribution of 
the cc resonances is judiciously subtracted from the latter on the experimental side. Since 
the NNLO formulae for the matrix elements given in [39J are valid only for s < 0.25, 
following the prescription of ref. [2H1 we have used for the region s > 0.25 only the 
formulae of ref. [SI] with = 2.5 GeV (because for s < 0.25 the formulae of [21] with 
fib = 2.5 GeV quite accurately reproduce the full NNLO results obtained with fib = 5 
GeV) and assigned to the integral over this range of s the same fib uncertainty as has 
dr{13 — i> XJ,^l~)/ds computed for s = 0.25. We get in this way 

BR{B ^ X,/i+/i")nonrcs = (4.39+°;^^ ± 0.24) x 10"^ , (16) 
BR{B ^ X,e+e-)nonres = (7.26;°:| ± 0.28) x 10"*^ , (17) 
^This conclusion has been reached in a discussion with M. Misiak. 

^Keeping TOs 7^ has numerically a very small impact on dT/ds itself but Smax < 1 for the upper 
integration limit partly cures the problem associated with the nonperturbative contributions to the dif- 
ferential rate, which for s 1 dominate in the expression and make it negative in the vicinity of 

s = 1 nn. 
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where the meaning of the errors is as previously. Taking into account the remaining 
uncertainties we estimate the total uncertainty of BR(B — > XsfJ.^ names for li4% and of 
BR{B -'fse^e~)nonres for -i4%- Our central values are in good agreement with the ones 
given in ref. [2Hj but due to the normalization to the width of the semileptonic charmless 
decay the overall uncertainty is smaller even though we take into account the uncertainties 
related to the electromagnetic correction. Within the errors and uncertainties the SM 
prediction ()16p is roughly in agreement with the published BELLE and recent BaBar 
results, which together give ^S] BR{B -^s^^^~) nonres = (6.2 ± 1.7) x 10~^, averaged 
over I = n,e, for the dilepton invariant mass i/g^ > 0-2 GeV. ^ 

For a relatively clean comparison of the experimental measurements with the theoret- 
ical predictions of interest can be also the rate integrated over the region of s above the 
cc resonances. We get there 

BR{B ^ X,/+/-)o.65<s<s._ = (2.32l[!i^ ± 0.14) x 10"^ (18) 

{I = e or /i) where the first uncertainty, corresponding to the fib dependence, is estimated 
with the help of the prescription of ref. [2H] described above. Better estimate of this 
uncertainty will become possible once the calculation of ref. PU'] is available. However for 
this range of s the nonperturbative I /ml corrections of refs. [37,.38j constitute yet another 
potential source of uncertainty. For s > 0.8 these corrections cannot be calculated reliably 
PH] {sm = 0.65 of that paper corresponds to s ~ 0.8) which manifests itself in the negative 
values of the expression © for s — 1. To estimate the uncertainty introduced by this 
factor we have computed BR{B -^s/^/^)o.65<s<smax switching off the l/ml corrections in 
© fors > 0.8. At/ife = 2.5 GeV this gives Si?(S ^ X,/+/-)o.65<.<wx = 2.66x10"^. The 
difference of order 15% between this result and ()18|) can be interpreted as the uncertainty 
associated with the 1/m^ corrections. Adding all uncertainties in quadratures we finally 
assign to the result (fTH|) the uncertainty of order ±20%. 



3 Scalar flavour changing neutral currents 

Even in the MFV MSSM with tan /3 ^ 1 ordinary one loop corrections involving charginos 
and stops can generate substantial FV couplings of neutral Higgs bosons to the down- 
type quarks {q = s,d) El 13 112] • For sparticles sufficiently heavier than the charged 
Higgs boson (which sets the mass scale of the MSSM Higgs sector, as in the MSSM for 
Mh+ ~ 200 GeV Mh ~ Ma ~ Mh+) the effects of these FV couplings can be described 
by the local Lagrangian of the form: 

= -QL [Xlr]"' 6ij(//° - 2A°) - qn [X^^Y' bL{H' + zA') + H.c. , (19) 

^Our result for BR{B Xse+e~)„o„rcs for ^/f > 0.2 GeV is similar to JTHl: (4.48i^:37±0.24) x lO^^. 
In the comparison with the BELLE result one has to take also into account the error in translating the 
"reduced" invariant mass s = q'^/ml into the experimental cut on physical q^. 
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where in the so-called approximation of unbroken SU{2) x U{1) symmetry the amplitudes 



[Xlr]''^ are given by 



mt, fmt\^ tan^ p V*Vtb 



^ ^ AMwKMwJ (l + ib tan (3) {1 + eo tan /3) ^' ^ ^ 

The factors ey ~ 0{1/16it^), cq and ib (see ref. [20] for the analytical expressions) depend 
on sparticle mass parameters; in particular, ey is directly proportional to the mixing of 
left and right stops, that is to the parameter At fl^. The factors eo and e?, which depend 
on both, as and the top Yukawa couplings, ensure proper resummation of the (tan/?)- 
enhanced terms from all orders of the perturbation expansion jHl ^1 EOl 1^ ■ Their 
signs and magnitudes depend directly on the signs of the supersymmetric /i and param- 
eters. Generally, the resummation factors suppress the FV couplings for fi > J3] and 
enhance them for fi < ^517. The amplitudes [Xrl]'''^ of the transitions 6^ —* Sfi{dfj) are 
given by similar expressions but with mi, replaced by rris^d) and are, therefore, suppressed 
(but are, nevertheless, important for the B^^-B^ mixing [HI El CH|)- The approximate 
formula (j^UI) captures the main qualitative features of the FV couplings generated in the 
MFV MSSM. For more accurate estimates of their magnitude and dependences on the 
MSSM parameters one has to use, however, more complicated approach developed in ref. 
|2()j which combines the resummation of the (tan/3)-enhanced terms with the complete 
diagramatic 1-loop calculation. In principle, for Msusy ^ My/ one should also take into 
account that the couplings (fTUj) are generated in the process of integrating out heavy 
sparticles at some scale /i^ ~ Msusy and should be evolved down to the matching scale 
/io using the RGBs similar to the RGEs for the quark Yukawa couplings in the SM 

where we have retained only the effects of the QCD renormalization. As a result, the 
couplings [Xlr]'''' would be multiplied by the factor [a.silJ'o) / Ois{iJis)Y^'^ , equal (for /xq = rrtt) 
1.073 for Us = 500 and 1.12 for = 1000 GeV. To take consistently such effects into 
account one would have also to determine sparticle couplings at the scale fis (and use them 
to compute the amplitudes [Xlr]''*). Since for the correlations discussed in section 4 only 
the values of [^lr]"^^ at /iq matter we will simply assume that sparticles are integrated 
out at the same scale /iq = mj. 

With the FV couplings (j^UI) the tree-level exchanges of and A^ generate at the 
scale /io Wilson coefficients of the and Op operators 

gl rrumf^ifio) / m, \3 tan' ^ V*Vtb ^ 

^^'^ - "8M1 [m^J (l + e,tan/5)(l+eotan/5)'^ ~ "^^ ' ^^^^ 

Note that the expressions for Cg^^ and Cp^^ through their dependence (via eo and e;,) on 
the coupling constants and at = vH^t^ (where yt is the top-quark Yukawa coupling) 
resumm terms of order tan"+'" (3 {n,m > 0) from all orders of perturbation theory. 
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Since the operators mi,Os,p, are renormalization scale invariant with respect to the strong 
interactions, the QCD evolution of 



r<i - W(i) , W(2) 



'S,P 



Air 



S,P 



+ ■ 



(23) 



reduces to the multiplication of Cg p{fio) by the factor [mb{fib)^^/m^^{fio)]. If as in 



Cgp(/io) oc m^^(/io) the dependence on of the formula for BR{B^ l^l ) cancels 
against the factor l/m^^(/i;,) present in the matrix element of the Os,p operators j4i 



Ml, 

q 



mr^(/i,)+mMS(^,) 



m 



MS 



(24) 



Complete 0{as) calculation of the scalar operators contribution to BR{B^ l^l ) in the 
MSSM would therefore require only computing higher order corrections to the matching 

1(2) l(2\ 

conditions at the scale /xq, that is to resumm all contributions to Cg and Cp of order 
Q;s(agaJ" tan""*"™ /3) for n, m > 0. 



One can also take a more general point of view and assume that the scalar operators 
O^g p are generated at the scale /xq by some yet unknown physics and investigate their 
effects on the h —>■ sl~^l~ and b —>■ dl^l~ transitions without any reference to the more 
fundamental theory, treating the Wilson coefficients &g p as free parameters. Assuming 
dominance of the scalar O^p operators, the formula for r(i?° ^21 E] takes the 

form 



r(RO ^ /+/-) ^ M«o 



647r3 



MS 



\v;M 



that is 



5i?(5° ^ l+r) ^ 4.27 X 10" 



1 

X- 

2 



The recent CDF upper limit 
therefore the stringent bound [IS] 



C'sif^b 
BR{B's - 



238 MeV 

2 



+ 



0.04 

2' 



4.2 GeVV 



MS/ 



(25) 



^1+^1-) < 0.95 X 10-6 at 90 % C.L. sets 



C's{^^b 



+ 



C'^(/i6 



< 2.2 X 




(26) 



Similar bound can be also derived for Cg p{^b) by using the corresponding experimental 
upper limit BR{B^ e~^e~) < 5.4 x 10"^ [45j but it is two orders of magnitude weaker. 
Analogous bounds on the (universal under the assumption of MFV) Wilson coefficients 
Cgp{fib) that can be derived from the experimental upper limits BR{B^ —>■ fi~^fj.~) < 
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1.6 X 10-^ [23 123 and BR{B^ e+e") < 8.3 x 10"^ 03] are less interesting as they 
depend on the value of \Vtd\, determination of which can be also affected by the new 
physics that gives rise to the scalar operators 



As follows from the formula ()22|1 . 



in the MSSM C's p 



■1(1) 



S,P 



oc mi, so that the 

effects of the scalar operators can be measurable only for the ii^ ^~ and t^t~ channels 
(the latter being very difficult for experimental searches so that at present no limit on 



T ' T 



is available). For tan/? ~ 40 — 50, substantial stop mixing and /i < 0, 



when the resummation of the leading tan" (3 terms enhances the FV violating couplings, 



IC: 



Ml) I 



|(-<Mi)| 



could be as large as ~ 10 leading to BR{Bl fi+fi-) ~ lO'^ ^U^. 
The bound eliminates therefore a large portion of the general MSSM parameter 
space. Moreover, as has been demonstrated in ^71 1201 ' such cases also the contribution 
of the FV couplings of the neutral MSSM Higgs bosons to the B^-B^ mass difference AMg 
is large and the experimental limit [AMgY^^ ^ 14/ps [26^ becomes in most cases more 
constraining (see the next section). 



It should be stressed, however, that the the bounds like ()26|) are completely indepen- 
dent of the specific way of generation of the coefficients jC'gpl and are valid generally, 
and not only in supersymmetry.^ In particular, one can imagine that the operators Ogp 
are not due to to the neutral Higgs boson exchanges between the FV violating down-type 
quark vertices and the leptonic vertices in which case sizeable effects of the scalar opera- 
tors Ogp could be present in any of the b — s{d)l^l^ transitions (for any lepton) and not 
accompanied by large contributions to the B^-B'^ mixing amplitude as in the MSSM. For 
this reason, the remaining analysis of this section will be done in a general framework. 
We will return to the MSSM only in the next section. 

The general bound ()26|) on jC^pl allows for an immediate estimate of the impact, 
the scalar operators Ogp may have on the rate of the inclusive process B — >■ XsH~^n~. 
Similar estimates can be also made for B X^e^e" and B XsT'^t' processes. From 
the formula (jH)) for the contribution of the scalar operators to the differential rate we get 
da 121]: 



d 

ds 



ABR{B 



3 



BR{B ^ X^eue) 1 



2qs (rrif,) 



2 3 



Cp{^b 



27r 



+ 6-^C^^(/i,)Cit(s,/ife) 



(27) 



where we have used the normalization to the width of the semileptonic charmless decays 

''The bound H26fl is valid also if the new physics, which gives rise to nonzero p involves sources 
of FV other than the CKM matrix, provided the coefficients C\ p are (superficially) normalized as in 
(jJJ. More generally, for a given lepton pair l^l^ the experimental upper limits on BR{B'^ l^l^) and 
BR(B^ set then independent bounds on the products (assuming that the Wilson coefficients 



are stiU normalized as in ll^t^Vftp 



&P 



for q ^ s and q = d, respectively, which can be 



directly used to constrain the maximal possible effects of the scalar operators in inclusive or exclusive 
B XJ+l- and B XJ+l' decays. 
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and for simplicity dropped the nonperturbative correction factor appearing in the denom- 
inators of the formula ()13|) . As remarked below the formulae (fT^ . the contribution of 
the operators to the inclusive rate BR{B Xs^~^fi~) depends on the choice of 

the renormalization scale fib- Since following ref. jH] we use ml^ = 4.69 GeV leading to 
m^^{m^^) ~ 4.2 GeV for the value of the running 6-quark mass, in what follows we will 
treat as free parameters the Wilson coefficients p taken at fit, = 4.2 GeV. The uncer- 
tainty related to the variation of the scale fi^ — > fi'b in the formula (j77|) is then roughly 
(ascribing for the estimation purpose to the interference term the same Hb dependence 
as have the other two terms) given by [mf^^{fi'h)/m^^{fib)f and is estimated to be 
This uncertainty has to be, of course, combined with the ones stemming from unknown 
electromagnetic corrections and the C-factor p4j) . Inserting numbers in the formula ()27|1 
we get 

^ABR{B X,/i+/i-) ^ 4.7 X 10~^ 
as 



C^ifib) ' + s C^pifib) ' + i'^C^p{fib)cS] . (28) 

Integrating over the full (0, 1) range of s and taking into account the limit ()2(jj) with 
m^^dib) = 4.2 GeV for /i^ = 4.2 GeV we obtain the estimate of the maximal possible 
contribution of the scalar operators to the "non-resonant" branching ratio: 

ABR{B ^ X,/X+/i~)nonres < 1.7 X / X (^1 ± 0.5 X X 10"' (29) 

where 

o.78</< 1.32 ,30) 

and the factor 

< r < 2 (31) 

depends on the relative magnitudes of \Cp\ and \Cg\: r = for \Cp\ = and r = 2 for 
\Cg\ = 0; for \Cp\ = \Cg\, as in the MSSM, r = 1. The ± refers to the two possible signs 
of the interference term depending on the sign of Cp (the interference is constructive for 
Cp < 0). We have used the approximate SM value C*^o(s, yUfc) ~ cIq' ~ —4.2 and mb = 
^poic _ ^ g ^Yie interference term. Thus, the maximal effect of the scalar operator 

is 3.7 X 10"^ for / = 1.32 and r = 2 (2.55 x 10"^ for / = r = 1). Comparing with the SM 
result |T6|l we conclude that the maximal contribution of the scalar operators allowed by 
the CDF limit (0) is at most at the level of 8% for / = 1.32, r = 2 (5% for / = r = 1), 
that is, substantially smaller the estimated uncertainty of the SM prediction. This is in 
sharp contrast with the findings of ref. where it has been claimed that even within 
the so-called minimal SUGRA framework the ratio BR{B^ fi~^ ^~)/ BR{Bg e^e~) 
can reach values as big as 2-3, corresponding to the contribution of the scalar operators 
as large as 100-200%. 
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For the branching ratio integrated over the range 0.05 < s < 0.25 we find 



ABR{B ^ X,/iV")o.o5<.<o.25 ^ 0.43 xfx[l± 0.88 x y^j x 10"^ (32) 

that is, the maximal effect is again of order 8% for / = 1.32, r = 2 (5.5% for / = r = 1), 
much smaller than the estimated uncertainty of the SM prediction for this range. For the 
range of s above the cc resonances the limit (PJ implies: 

ABR{B ^ X,/x+/i")o.65<.<i ^ 0.22 x / x (^1 ± 0.16 x x 10^^ (33) 

For this s range the maximal possible contribution of the scalar operators increases the 
branching fraction by ~ 15% for / = 1.32, r = 2 (11% for / = r = 1), that is again the 
effects of the scalar operators are not greater than the estimated uncertainty of the SM 
prediction.* Estimates of ABR{B Xse~^e~) can be also obtained in a similar manner. 

Experimentally first measured were the exclusive B decay modes B Kl^l~ and 
B K*l^l~ jH]. For B Kl~^l~ , which will be of interest for us here,^ the recent results 
for the "nonresonant" rates are BR{B Kfi+jj,-) = {A.StH ± 0.3 ± 0.1) x 10"^ 

and BR{B KI+1-) = (A.StH ± 0-3 ± 0.1) x 10"^ averaged over e and /i (BELLE) and 
BR{B Kfi+fi-) = (A.Stli ± 0.4) x 10"^ and BR{B Kl+l') = (6.9ll| ± 0.6) x 
10^^ (BaBar). The main uncertainty of the theoretical BR{B Kl^l~) calculation 
is related to the determination of the nonperturbative matrix elements of the relevant 
operators between the initial and final meson states. Different techniques used for this 
purpose resulted in the SM predictions for this branching fraction spanning the range 
(3.0 - 6.9) X 10"^ HZl UHl EB]. Within the experimental errors the new experimental 
results are in fair agreement with the SM-based NNLO theoretical estimate given by Ali 
et al. BR{B Kl+l-)nonres = (3.5 ± 1.2) x 10"^. Substantial lowering of the SM 

prediction compared to the earlier one of Ali et al. (based on the NLO calculation) 
BR{B Kl^l~)aonres = (5.7 ± 1.2) X 10^^ was mainly due to the superficial lowering of 
values of the formfactors parametrizing the operator matrix elements. This was motivated 
by the fact that the = value of the Ti(g^) formfactor obtained using the so-called 
QCD light cone sum rules (LCSR) gave, compared to the data, too high a branching 
fraction for the B — > K*'y mode jlH], suggesting that the LCSR method systematically 
overestimates the formfactors. 

The contribution of the scalar operators to the branching fraction BR{B Kfi^fi^)noni 
has been analyzed in ref. ^H]- At that time only the upper limit BR{B^ K^fi^fi^) < 
5.2 X 10"^ was available [2Z], so the conclusion of ref. was that the constraint imposed 
on l^^p + \C^\^ by the limit BR{B'^ -> /i+yU~) < 2.6 x 10"^ was significantly stronger 
than the one that could be obtained from the limit on BR{B — K~^fi~^fi~). With the 
new numbers the situation is somewhat different and we summarize it below. 

*With the old limit BT{B^ — > fJ,^^J.^) < 2 x 10~^ the effects of the scalar operators in this range 
of s could be almost twice as big as the estimated uncertainty. 

^As analyzed in ref. the contribution of the scalar operators to BR{B K* ji^ ji^) is too small 
to be interesting. 
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The scalar operators contribution to the nonresonant branching ratio can be written 
as [131 



ABr(B Ktl 



cm 



X {q^fi][q^)\bFs\^ + q^\bFp\^ + 2g^Re(F;5Fp) 
+2mz(M| - M| + g2)Re(F;5Fp)} 



(34) 



where q^ is the physical lepton pair invariant mass, /^/(g^ 



1 — 4m^ / g2 and 



5F. 



5,P 



iap(/i,) 



(/^fe) 



MS 



(M| - Ml) /o(g^ 



(35) 



miC 



off 
10 



Af| - Ml r 




The coefficient differs from C^o (s, /i^) given in eq. (fT^ by setting to zero the functions 
Lt)o,{s) (the effects of (-019(5) are supposed to be taken into account in the formfactors /o(q'^) 
and /+(q'^)). Note that CJ^g i^J; hence the whole formula ((311), independent of 
the renormalization scale /Xb. Following the recipe of ref. |2H1 for the central values of 
the formfactors foiq"^) and f^{q^), as well as for fri'f') appearing below, in eq. ((33), 
we use their lowest values obtained within the LCSR approach which amounts to using 
the formula (3.7) of jlB] with the parameters collected in Table V of that paper. At 
the same time, again following ref. we ascribe to the values of the formfactors 

the uncertainty of order 15%. The formfactors introduce therefore in the results for 
{d/dq'^)ABr{B -^^^/^)nonres the largest (barring the discussion how big errors are 
introduced by using the effective Lagrangian with non-local coefficients Cg^(g^), Cj^{q'^), 
for the exclusive process) uncertainty of order 30%. 



Integrating over in the kinematical limits 4m^ < q^ < {Mb — MkY ^ind assuming 
that the new physics contribution to Wilson coefficients other than C^p is negligible we 
obtain for the dimuon mode 



6.36 X 10 



-8 



X <^ a 



+ 



(36) 



where for = 4.2 GeV a = 0.30 ± 0.10 and b = 0.19 ± 0.06. The uncertainties of a and b 
are due to the uncertainties of the formfactors /o(g^) and f+{q^). Through the formfactors 
the total uncertainty of the scalar operators contribution is obviously strongly correlated 
with the uncertainty of the SM prediction. Sticking to the central values of a and b and 
taking maximal values of \Cs\ ^'^d \Cp\ allowed by the bound (^^1 we get 



< 0.8 X / X 1 ±0.45Jr// x 10 



1-7 



(37) 
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that is, the maximal possible contribution of the scalar operators to the nonresonant 
branching fraction can be (for Cp < 0, C*^ = 0, and the lowest possible value of Fb^, i.e. 
for the + sign, r = 2, / = 1.32) as large as 1.7 x 10~^, roughly of the same magnitude as 
the error of the experimental result and 1.5 times bigger than the estimated |48 [ I13 | l28]) 
uncertainty (~ 1.2 x 10"'') of the SM prediction. Similar estimates can be also done for 

ASr(S ^ Ke+e-)nonres. 

Finally, an experimentally interesting quantity |13| I3U] may be the integrated over 
forward-backward lepton asymmetry measured in this decay given by^° 



™ ~ BR{B Kfi+fi') V 167r2 ) 7rM| 



2 



X J dq^ miX{q\ Ml, Ml) Pf{q^) ReiF*SFs) (38) 



where 



Fv = Cf{^^,)U{q') + 2m,Cf(/.,) ^^^|^2f^ (39) 

with Cg^ and C^^ differing from Cf^ and C^^ of eqs. ()12|1 by setting to zero^^ the functions 
ujg{s) and 0^7(5). The asymmetry ApB vanishes in the SM in which Fs = SFs = 0. For 
the dimuon channel, integrating over the whole range and using Hh = 4.2 GeV we get 



< ± 



4.8 X 10 



-7 



BR{B K/j+fi- 



X (1.5 ± 0.5) X Jr'f % (40) 



where < r' < 2 (r' = for = and r' = 2 for = 0; for \C^\ = \C^\, as 
in the MSSM with, r' = 1). The uncertainty of this result being dominated by the 
30% uncertainty arising from the formfactors f+{q^) and fri.Q'^)), is of course strongly 
correlated with the uncertainty of the total branching ratio. Still, the maximal possible 
asymmetry allowed by the limit (j26|) is of the order of a percent and may be detectable 
in the future. 

We conclude that given the experimental limit IQ, the effects of the scalar operators 
in the inclusive process are typically of order 5 — 15%, always smaller than the estimated 
uncertainty of the SM NNLO prediction. On the other hand, the maximal allowed contri- 
bution of the scalar operators to BR{B — > Kfi~^^~), although larger than the estimates 
of the theoretical uncertainty of the SM prediction made in jlHl UHl I2H] , is only roughly 
of the order of the present experimental error. While the latter can shrink in the near 
future, the spread of the different SM based theoretical predictions and the problems 
with the formfactor values obtained using the QCD LCSR may suggest that the true 

-'^^ Of interest can be also unintegrated differential asymmetry 
"'^^for q^/mf > 0.25 we also set to zero the functions fI'^'^\ 
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uncertainty of the SM prediction is larger than estimated in [JHl ^1 12H] , thus preventing 
the rehable comparison of the theoretical predictions with the data. The forward back- 
ward asymmetry of the muon distribution, if detected in the high statistic data, could be 
also indicative of the scalar operators contribution (the asymmetry vanishes if only the 
SM operators contribute) but its translation into the values of Cg and Cp depends on 
the formfactors too. Thus, before the status of the formfactors is clarified and the er- 
rors associated with them reliably estimated the exclusive mode B K^'^n~, although 
potentially interesting, will not be able to put constraints on the coefficients and Cp. 

The coefficients and Cp of the most interesting (largest allowed) magnitude cannot 
be however, as in supersymmetry, due to the tree level exchanges of the neutral Higgs 
bosons between the effective quark FV vertices and the Higgs-lepton-lepton vertices. As 
we shall see on the MSSM example in the next section, in such a case possible effects of 
the scalar operators (apart from being slightly reduced by the relation jC^^'*! ~ |C*p^''| so 
that r = r' = 1) can be further constrained by the B'^-Bg mixing. 

4 Correlation with BR[B Xs^) and the B^-B'^ mass 
difference 

In assessing potential effects of the scalar operators in the preceding section we have 
ignored the fact that the new physics, which gives rise to them, can also modify the 
remaining Wilson coefficients. In the MFV MSSM charginos and stops as well as the 
charged Higgs boson contribute to Cio(yUo) and Cg{fio) through the box, Z°-penguin 
and, in the case of Cg{fiQ), also through the photonic penguin diagrams. Likewise the 
coefficients of the Cj and Cg are modified by loops containing these particles. It should 
be also stressed that supersymmetric contributions to Cj^ (/io) in eq. (jlj) necessary for 
complete NNLO calculations are only partly known for C7 and Cs (and only for a sce- 
nario with light right-handed stop and charginos) jS^j and are unknown for the other 
coefficients in (j2)). Out of the relevant for the b — > sl~^l~ transition Wilson coefficients 
only the modulus of C!}^ (but not its sign) is rather well constrained by the measurement 
of BR{B — > Xsj). The other coefficients can still accommodate substantial new physics 
contributions. 

If is light - a necessary condition for generating in the MSSM nonegligible Wilson 
coefficients of the scalar operators - its contribution to C^^ is substantial and has the 
same sign as the SM contribution. Therefore it must be cancelled out by the chargino- 
stop contribution. For tan/3 ^ 1 the latter is proportional to tan/3 and can be very 
large if these particles are light. Its sign depends on the sign of Atfi and for A^fi > (in 
our phase convention) it is opposite to the sign of the W~t and H^t loops so that the 
cancellation is indeed possible. Since for At^ > the Wilson coefficient Cp^^ is negative 
{Cg^^ is positive), the requirement that the calculated in the MSSM BR{B —>■ Xg^y) agrees 
for a light with the experimental result necessarily leads to positive contribution of 
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the interference term Re(C5'Q Cp^"**) in the formulae Q and (fTfj) (recall that the SM 
contribution to C[q' is also negative) so that in the estimates (j29|) . (j32|l . (j33|l . pTjl and 
(jinjl the + signs apply. 

In principle the chargino-stop contribution could even reverse the sign of C^^ leading 
to a value of BR{B Xg'j) compatible with the experimental result. The sign of the 
Ke{C!}^CQ^*) term in the formula IQ would be then changed modifying predictions for 
the B — > Xgl^l^ rates. Such situation, which could most easily be distinguished by mea- 
suring the dilepton invariant mass spectrum and the forward backward asymmetry in the 
BR{B — > Xsl~^l~) [Sn], requires light, ~ 100 GeV, charginos and stops and, for light 
and tan /9 3> 1, is strongly fine tuned 13^. Much more natural appears the possibility that 



charginos and stops are rather heavy and their contribution to Cj^ , despite substantial 
stop mixing necessary for generating large Cgp, is small, just of the right magnitude (and 
sign) to cancel the contribution of the charged Higgs boson. In such a scenario the value 
of C^^ must be close to the one predicted in the SM and the contributions of stops and 
charginos to C{o''(yUo) and Cg^\fiQ) is, as we have checked by using the formulae of ref. 
negligible. 

The contribution to Cl^'lfio) and Cg^\no) through the box diagrams, and 
photonic penguins has been computed in ref. [33]. For tan/5 ^ 1 these contributions 
are not enhanced and are negligible even for the charged Higgs boson mass as low as 200 
GeV. As has been found in [56^ 12} the H^t loops also generate the FV couplings (fT^ 
and the resulting contribution to p grows as tan^ (3. However, for Mh+ ~ 200 GeV and 

tan /3 ^ 50 this contribution to the coefficients Cg^p are roughly two orders of magnitude 
below the upper limit (pUj) and, hence, their impact on the B Xsl^l~ rate can also be 
neglected. 

Thus, for sparticles heavier than, say, 500 GeV, the only sizeable SUSY effects in the 
h — > sfi~^fi~ transitions can be due to scalar operators. 

As has been observed in [^l El 1201; MFV MSSM whenever the coupling [Xlr]'''' 

()20|) is large, the tree level exchanges of the neutral Higgs bosons if ° and between the 
tree-level effective vertices fll9p give also large negative contribution to the mass difference 
AMg of the -Bg and B^ mesons. In the so-called approximation of unbroken SU{2) x U{1) 
symmetry, in which also the formula (j^n|) is valid, one gets [20j 



12.8 
ps 



tan P 



50 



238 MeV 



ts 



0.04 



mfe(/io 



3 GeV 



60 MeV 

1 2 



'1 + eo tan l3) (1 + eb tan (3) 



(41) 



(the analogous contribution to the B^-B^ mass difference, being suppressed by the ratio 
^d/^si is negligible). Typically the couplings [Xlr]^^ which give rise to Cg^p saturating 
the bound (PU)) lead to AM^ below the present lower experimental limit ~ 14/ps 
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Figure 1: Scatter plots of (1/2) (|C^^^^p + iC^^^^p) versus AM, in the MFV MSSM 
for sparticle masses greater than 500 GcV. Panels a-f correspond to {Ma, tan f3) values 
(200,40), (200,50), (300,40), (300,50), (400,50), (500,50), respectively. Points to the left 
and above the solid lines are for Fb^ = 238 MeV excluded by AM, > 14/ps and BR[B^ — > 
/X+/X") < 0.95 X 10~^, respectively. The same constraints for Fb, = 207 MeV and 269 
MeV are shown by dashed and dotted hues, respectively. 
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In order to see how the possible effects of the scalar operators in the b — >■ s/x^/i tran- 
sitions are limited by the experimental lower bound on the B^-B^ mass difference AMg we 
present in figs. [T^-UlF scatter plots of the combination (1/2) {\Cs^^\iJb)\'^ + |Cp ~ 

\Cg^^\fib)\'^ ~ \Cp^^\iJb)\'^ for jjh = 4.2 GeV versus AMg calculated using the approach 
developed in ^J] for a few combinations of the parameters {Ma, tan (3). The plots 
have been obtained by scanning over the MFV MSSM parameters (in the sense explained 
in more detail in ref. [201) with the lower bound on sparticle masses Msusy ~ 500 GeV. 
More specifically, we have scanned the relevant parameters in the ranges such that: 500 
GeV< mc, < 1 TeV, with 0.75 < |M2//i| < 1.5, 1 < m^/Ma < 3; 0.7 < Mijmc^ < 1.3, 
1.1 < MijMi^ < 1.7 and -35" < 0^- < 35"; 0.5 < Mi^/m-g < 0.9, with A = A; masses 
of the first two generations have been taken as max(M^^, M^^). All points, for which 
computed BR{B Xg'j) does not agree with the experimental result have been rejected. 
We have used \Vts\ = 0.04 and Fbs = 238 MeV but the limits for other values of these 
parameters can be obtained by simple rescalings. Horizontal lines in figs. show the 

upper bound ^ on (1/2) {\C^\^ + jC^p) for /if, = 4.2 GeV and Fb, = 238 MeV (sohd 

lines), Fb, = 207 MeV (dashed lines) and Fb, = 269 MeV (dotted lines). Vertical lines 
show the corresponding constraint imposed by the experimental lower limit AMg > 14/ps 
(excluded are the points to the left of these lines). 

From figs. ^-QF it is clear that the lowest possible values of Fb^ which in the 
model independent analysis of the preceding section gave the biggest effects of the scalar 
operators in the B — > Xg^^ ^~ and B Kfi'^^~ transitions, are in the MSSM allowed 
only for small values of the Wilson coefficients IC^^^"*!^ and ICp*"^"*!^. Moreover, for Ma > 
200 GeV and Fb^ ~ 238 MeV the lower limit AMg > 14/ps becomes more constraining 
than the bound ([Q). This means that in the MSSM (or any other model in which Os,p 
arise from the FV couplings similar to [Xlr]"^* in eq. (HHI)) the possible effects of the 
scalar operators Os,p in the B XsfJ^^fJ^^ and B Kfi^^~ decays must be smaller 
than the estimates given in section 3. For example, using the formulae of section 3 and 
the numbers that can be extracted from fig. ^p, we find that for M^ = 300 GeV and 
tan 13 = 50 the maximal effects in the inclusive process are bounded by 

ABRiB ^ X,/i+/i-)„s ^ 2.2 X 10"^ (42) 

obtained for / ^ 12/14 (for which and AMg > 14/ps allow for maximal value of 
\Cg\ ~ \Cp\) and r ^ 1, that is, any effects of the scalar operators must be below 5%. 
The maximal effects in the exclusive decay B — > Kfi'^fi~ are then also suppressed by the 
limit on the B^-B^ mass difference: 

ABR{B ^ Kfi+ fi')^onres ^ 1-0 X 10"' . (43) 

The suppression further with decreasing value of tan/3 and increasing mass scale of the 
Higgs boson sector (set by Ma) up to Ma ~ 650 GeV. 

^^In producing these plots we have corrected a bug in our fortran code which resulted in using in 
refs. [nilTHlEni C^%\fJ'o) instead of C^^p\tib) in calculating BR{B° ^ A^+M")- As a resuh numer- 
ical values of this ratios in figures of these references should be rescaled upwards roughly by a factor 
[mp(4.2 GeV)/mfS(TOj]2 ^ 2.36. 
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Since the effects of tlie FV couplings (fT^ in BR{B^ fi~^fi ) scale as (l/M^)^ while 
in AMg only as (l/M^)^, for sufficiently heavy Higgs sector and sufficiently large couplings 
[Xlr]''' it is possible to get from the formula (gH) 6{AM,) < 2{AMsf^ that is |AM,|^ssm 
again compatible with the experimental lower limit (this possibility is seen in the upper 
branch of points in figure Qf) and, at the same time, BR{B^ A^^A* ) below the CDF 
upper limit. This happens only for Ma ~ 750 GeV. For such Higgs boson masses and 
values of the couplings [Xlr]'*'' the upper bound can be saturated and simultaneously 
Fb^ can assume lowest possible values obtained from lattice simulations Only then 
could the effects of the scalar operators Os,p in B ^ X^yU+yU^ and B — > Kfi^^^ decays 
reach the maximal values discussed in section 3 (reduced only slightly by the fact that 
in the MSSM r = r' = 1). One should stress, however, that, at least in the MFV 
supersymmetry, the couplings [Xlr]**'' of the required magnitude can be generated by the 
chargino stop loops only for very large values of the stop mixing parameter At along with 
significantly split stop masses and are very unlikely from the point of view of generation 
the soft SUSY breaking terms and most likely leading to the dangerous color breaking 
minima of the scalar fields potential. 



5 Conclusions 



Rare decays of B mesons are one of the places where the ongoing experimental mea- 
surements can reveal effects of new physics. The processes involving the b sl^l~ and 
b dl^l^ transitions are particularly interesting in this context. The most general low 
energy Hamiltonian describing their phenomenology involves the so-called scalar opera- 
tors (9g = {sLbR)(Jl), Op = {sLbji)(I'y^l) (and similar ones with sl ^ cIl). Their Wilson 
coefficients are negligible in the SM but in models of new physics can be quite substan- 
tial compared to the coefficients of the other, usually studied, operators. This is so for 
example in the minimal supersymmetric extension of the SM even if the supersymmetric 
partners of the known particles are rather heavy, provided the ratio of the vacuum expec- 
tation values Vu/vd = tan/3 of the two Higgs doublets is large and the mass scale of the 
extended Higgs sector is not too high. 

In section 3 of this paper, following the earlier work J^]; we have used the experimental 
upper limits on the branching fractions BR{Bg^ to place the constraints on the 

Wilson coefficients of the scalar operators relevant for the b sl^l~ and b dl^l~ 
transitions. Particularly stringent constraint obtained from the limit BR{B^ f^^f^ ) < 
0.95 X 10"^ has been subsequently used to asses in a model independent way the impact 
the scalar operators {sibRj^fifi), (sibRj^fi'y^fi) may have on the rates of the inclusive 
B — > and exclusive B Kfi~^fi~ decays. 

We have found that the increase of BR{B — > X^yU+yU") due to the scalar operators 
cannot exceed (5-15)% (depending on the range of the dimuon invariant mass), that is, 
it is always smaller than the uncertainty of SM NNLO result which we have estimated 
in section 2. The large effects of the scalar operators found in this decay in ref. j^Hl 
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are therefore already excluded. On the other hand, the maximal increase of exclusive 
rate BR{B Kfi~^n~) can be still quite large, of order 1.7 x 10~^, comparable with the 
present error of the experimental result. The latter, when compared to the SM prediction 
BR{B Kfi~^fi~) = (3.5 ± 1.2) x 10^^ j2H], leaves some room for positive new physics 
contribution. However the SM prediction for this rate hinges on the theoretical problems 
related to the determination of the relevant nonperturbative formfactors. Before this issue 
is settled (and the experimental errors shrink) no firm conclusion about the detectability 
of new physics effects in the exclusive decay B Kn~^^~ can be drawn. 

In the supersymmetric scenario with large tan/3 and not too heavy Higgs sector, 
in which large values of the Wilson coefficients of the scalar operators can be naturally 
generated, the potential effects of Og and Op inb ^ sfi~^^~ are further constrained by the 
experimental lower limit on the B^-B^ mass difference. This has been illustrated in section 
4 in the case of the minimal flavour violation scenario considered in papers jHl El 
However, the limits on the scalar operator contributions to BR{B^g fi~^fi~), BR{B 
Xs/i^/i~) and BR{B —>■ Kn~^fi~) that can be derived by inserting in the formulae of 
section 3 numbers exctracted from figure ^ (for different values of tan/3 and M^) are 
valid also if the flavour violation originates in the squark sector, provided supersymmetric 
particles are heavy enough in order not to contribute appreciably to the box and vector 
boson penguin amplitudes. This is because the speciflc relation between the Wilson 
coefficients of Og and Op and the Wilson coefficients of the scalar operators contributing 
to the Bg-B^ mixing amplitude relies only on the existence in the low energy effective 
theory of the flavour violating couplings ()19|) and not on the speciflc mechanism of the 
flavour violation in the underlying theory. 

Note added While completing this paper we have learned about a similar independent 
study by F. Kriiger et al. jHT]- In particular they conflrm our conclusion that the large 
effects of the scalar operators found in ref. j^H] in the inclusive rate are already excluded 
by the experimental data. 
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